Recovery of Hidden Interference in Mott Insulators 
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Particle statistics plays a crucial role in a strongly interacting quantum many-body system. Here, 
we study the Hubbard model for distinguishable particles at unit filling. Starting from the superfluid- 
like state in the strong tunneling limit and gradually reducing the tunneling so that the on-site 
repulsive interaction dominates, the state ends up in a symmetric superposition of Mott insulator 
states. This result can be experimentally confirmed by the recovery of interference patterns in 
the density correlation functions. We also show that this state is a maximally entangled state, in 
contrast to the standard picture. 



The Bose-Hubbard model is a well studied model of 
spinless interacting bosons In this model, the bosons 
interact with on-site repulsive interactions and can tunnel 
between adjacent lattice sites. A quantum phase transi- 
tion occurs between a superfluid state in the strong tun- 
neling limit and a Mott insulator state in strong inter- 
action limit, as a result of the competition between the 
tunneling and the interaction [2]. Recently, this model 
has been experimentally realized using ultra-cold atoms 
in optical lattices and a phase transition has been ob- 
served by examining the diffraction patterns of the atoms 

i. 

The Mott insulator state was experimentally verified 
by the detection of a single Gaussian distribution in the 
diffraction pattern. Particles in this state are seemingly 
localized in each lattice site with vanishing spatial coher- 
ence over other lattice sites. In a more recent experiment, 
the measured spatial correlation function reveals higher- 
order coherence in this state [1, However, because 
of the symmetrization postulate for identical bosons, it 
is impossible to determine whether the observed inter- 
ference is a result of the trivial coherence due to the 
symmetrization postulate or is a result of the intrinsic 
coherence generated by particle-particle interaction. 

In this paper, we study a Hubbard model for distin- 
guishable particles (as opposed to identical bosons or 
fermions) to clarify the hidden coherence and entangle- 
ment in the Mott insulator state. We will show that 
quantum interference can be found for distinguishable 
particles, which demonstrates the (truly) coherent na- 
ture of the Mott insulator state, in which each particle 
not only occupies all sites in a linear superposition state 
but also has correlations with other particles. The theo- 
retical prediction of the recovery of the interference pat- 
terns can be experimentally tested by a relatively small 
scale experimental system. Using exact diagonalization 
methods, we derive ground state properties of the Hub- 
bard model for distinguishable particles. We show that 
the ground state of this model is a fully symmetric state, 
in complete analogy to its bosonic counterpart. The dis- 



tinguishability causes a large degeneracy in the Mott in- 
sulator states in the limit of zero tunneling. At small but 
finite tunneling, the ground state can be approximated as 
a symmetric superposition of such degenerate Mott insu- 
lator states. Quantum interference 0, appears in the 
correlation function for coincidence detection of particles 
because of this (non-trivial) superposition. A quantum 
many-body system with distinguishable particles can also 
be very useful in studying entanglement in the quantum 
phase transition, as it avoids the subtle issue associated 
with entanglement and quantum indistinguishability. We 
compare the entanglement in the distinguishable particle 
Hubbard model and that in the standard Bose-Hubbard 
model. We adopt an operational definition of the en- 
tanglement following Ref. [8], where entanglement is the 
correlation between modes (not particles) in subsystems. 
Very different results are obtained for the two models. 

Consider a Hubbard model for distinguishable particles 
in a one-dimensional periodic lattice. The particles can 
tunnel between nearest neighbor sites with tunneling ma- 
trix element to which is the same for all species of distin- 
guishable particle (henceforth called an "isotope"). Par- 
ticles located on the same lattice site interact with each 
other with a repulsive interaction Uq for all isotopes. In 
order to directly compare with the Bose-Hubbard model, 
we define an operator-based Hamiltonian for the distin- 
guishable model. The Hamiltonian can be written as: 

H = -toY^ + h.c}j + Y ^ n,(n, - 1) (1) 

a,i i 

where rii = Uai is the total number of particles on site 
and a labels the isotopes. The operators Ca,i and c]^ ■ 
are the annihilation and creation operators for isotope a, 
and satisfy the commutation relations [ca^i^ cj^ j] = ^cx(3^ij 
and [cq,^^, C/3j] = 0. When the index a is omitted, the 
model becomes the standard Bose-Hubbard model. 

We are interested in the case of unit filling for a finite 
lattice with a total of N sites, where the number of par- 
ticles equals the number of lattice sites. We also only 
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FIG. 1: Ground state properties with = 4. (a) Schematic 
energy structure of the ground states and the first excited 
states. Energy levels are shown as solid lines with its energy 
labeled on the left hand side and its degeneracy labeled on the 
right hand side. Energy spectrum versus to/Uo are plotted for 
the lowest 15 states in (b) the Bose-Hubbard model and (c) 
the Hubbard model for distinguishable particles. 



consider the case where there is one particle per isotope 
= Xli i^oc^i ~ ^' At to ^ ^0, whcrc the on-site in- 
teraction is negligible, each particle occupies the lowest 
single particle mode of the tunneling Hamiltonian. The 
ground state is a non-degenerate state with wave function 
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where cj^ = Xln ^*^^^a,n/V^ is the creation operator 
of the a single particle mode of momentum k for isotope 
a. This state, with an energy of —2Nto, corresponds 
to the superfluid limit of the Bose-Hubbard model and 
is symmetric under exchange of any two particles. In 
the first excited state, one particle occupies the second 
single particle mode. In an A/" = 4 lattice, the excited 
states are 8-fold degenerate with an energy 2to above the 
ground state as shown in Fig.[T](a). At to = 0, the parti- 
cles avoid each other because of the repulsive interaction 
and localize at different lattice sites. The ground states 
are degenerate states with A^!-fold degeneracy because of 
the distinguishability of the particles. The wave func- 
tions have the form = Yla=i ^i,i=p{a)\^) where 
the P{a) is a permutation of the list {1, • • • , N}. These 
states correspond to the Mott insulator state of the Bose- 
Hubbard model. With A/" = 4, there are 24 degenerate 
states, which are shown in Fig. [T] (a) together with the 
excited states. 

Exact matrix diagonalization is used to study the 
ground state properties of the above model. We calcu- 
late the eigenstate spectrum, the eigenfunctions, and the 
fringe visibility for a range of to/Uo. The energy spec- 
trum with = 4 particles is plotted in Fig. [1] (c). We 



find that the ground state energy and wave functions of 
the distinguishable particles are in exact agreement with 
that of the bosons (Fig. [11(b)) everywhere in the param- 
eter regime, except at zero tunneling. At zero tunneling, 
the statistics of the particles causes high degeneracy in 
the ground states. At to 7^ but to <C /7o, the ground 
state can be approximated as the symmetric superposi- 
tion of the Mott insulator states: 
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The symmetry in the Mott insulator state at to <C Uo 
is hence not a trivial consequence of the symmetrization 
postulate. Even though the gap between this ground 
state and the excited states is ~ t^/Uo^ very small at 
this limit, it can be prepared experimentally by an adi- 
abatic approach. Starting from the superfluid-like state 
at strong tunneling limit and slowly turning off the tun- 
neling to to <C /7o, the final state becomes the symmetric 
linear superposition state before it is thermalized by the 
environment. 

The visibility of the fringes in the single particle diffrac- 
tion pattern can be viewed as an order parameter for a 
finite lattice model and is defined as 



V 



max{(nfc)} - min{(nfe)} 
max{(n/e)} + min{(n/e)} 



(4) 



with Uk = k^oc.k being the number operator in k- 

space and k = 2i\n/N for n = 0, • • • , A/" — 1. It is an im- 
portant parameter for studying quantum phase transition 
[9]. As discussed above, at to ^ t/o, max{(n/e)} = A^ at 
A: = and min {(n/e)} = for all other /c, we have V = 1; 
while at to <C t/o, (n/c) = 1 for all /c, we have F = 0. Be- 
tween these two limiting cases, the visibility fringe mono- 
tonically increases from to 1, clearly demonstrating the 
change in nature of the ground state between the Mott 
insulator to the superfluid-like states (see also Fig. [2]) 
for the distinguishable particles. 

The detection of the correlation functions is a more 
powerful tool in studying the nature of a quantum many- 
body system than single particle detection 0]- In 
this process, the particles in their ground state are re- 
leased from the lattice and the spatial density correla- 
tion is measured after a sufficiently long time. To study 
this process quantitatively, we treat each lattice site as 
an harmonic potential, but with an energy barrier be- 
tween adjacent sites that is half of the harmonic poten- 
tial: Vo = mcj^ag/lG, where ao is the size of the unit 
cell in the lattice, m is the mass of the isotopes, and 
is the frequency of the harmonic potential. We describe 
the wave function of the localized state at site i as 

(x|cl^,|0) = </.(x-i?i) (5) 

with (\){x) = ^ J— e~^^^^^ in the Gaussian approx- 
imation. The position uncertainty of this state is 
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FIG. 2: Particle detection (single particle detection probabil- 
ity) versus position Xb at time t for the Mott insulator state 
(thick solid curve) and for the superfluid-like state (thin solid 
hne) with = 4. 

a = h/muo. Estimating the tunnelling amplitude us- 
ing to ~ huje~^^^^^ [10], the tunneling matrix element 
can be expressed as to = {Ji^ /ma)e~^o/'^^^ in terms 
of (J. Below, we assume Sh^/ma^ = 5/7o in the cal- 
culation. After being suddenly released from the lat- 
tice, the particles evolve freely with the transformation 

2 

U{t) = exp {—i ^^^) at time t. Here, Pa is the mo- 
mentum operator of particle a. In the coordinate basis, 
we have 

MUit)cU0) = (6) 

for particle a. 

The density operator at position can be written as 
n{xi)) = \xa = X})){xa = ^6 1 , including all isotopes. 
Using the wave functions derived above, the density of 
the particles at time t, also called the first order cor- 
relation, can be calculated as {t)n{xb)U{t)\ilj). In 
Fig. [21 we plot the density for both the Mott insulator 
state at to <C Uo and the superfluid-like state at to ^ Uo. 
For the Mott insulator state, the density shows a smooth 
distribution with a Gaussian profile resulting from the 
initial Gaussian state in Eq. ([5]). Note that an incoher- 
ent product state of localized states on all lattice sites 
has the same density distribution as plotted in the figure. 
For the superfluid-like state, an interference pattern ap- 
pears with three major peaks corresponding to particles 
traveling with wave vectors k = 0,±27r/ao respectively. 
This result agrees with the analysis of the visibility fringe 
defined in Eq. ^ for large N. 

Higher order correlation functions for the coincidence 
detection of particles can be calculated similarly. We plot 
the second order correlation function 

C2{xb) = {ij\U\t)n{0)n{xb)U{tm (7) 

in Fig. [3] (a) and the fourth order correlation function 

C^{xb) = {i^\UHt)n{0)n{x,)n{0)n{0)U{t)\i^) (8) 

in Fig. [3] (b). In the plots, the correlation functions 
for each (fictitious) Mott insulator state show the same 
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FIG. 3: (a) Interference in the normalized second order cor- 
relation function and (b) the normalized fourth order corre- 
lation function (lower plot) with N — 4 particles. Dashed 
curves: individual Mott insulator state; thick solid curves: 
symmetric superposition of the Mott insulator states in 
Eq. (|3)); and thin solid curves: superfluid-like state. 

Gaussian profile as their density distribution. The cor- 
relation functions for the superfluid-like state show the 
same interference pattern as their density distribution as 
expected (each particle is independent and uncorrelated). 
In sharp contrast, interference patterns can now be seen 
in the higher order correlation functions for the symmet- 
ric superposition of the Mott insulator state in Eq. ([3]). 
These correlation functions include linear superposition 
of terms given by Eq. (|6]) for all the permutations P. 
Quantum interference appears as a result of the coherent 
superposition and particle-particle correlation. It can be 
shown that for N particles, the interference pattern can 
be found from the second to the 7V-th order correlation 
function. 

In order to understand the above result, let us con- 
sider the simplest case of = 2. The Mott insulator 
state in this case is |?/;i2) = -^{\R)i\L)2 + \L)i\R)2). 

With |±)i,2 = ^(|L) ± |i?))i,2, \^12) = ^(|+)l|+)2 - 

| — )2)- If we ignore particle 2, the state of particle 1 
is a mixed state, pi = ^{\L)i{L\ \R)i{R\), and there is 
no interference. However, if we project particle 2 to |+)2 
by detecting it, the resulting state of particle 1 remains 
at -^(|L) + and there is interference. This is a 

standard situation of a "quantum eraser". This result 
clarifies the origin of the interference pattern (and hence 
coherence) in the Mott insulator state. As we are deal- 
ing with distinguishable particles, it is clear here that the 
interference originates from the particle interactions, not 
due to the symmetrization postulate. 

The above results are for small number of particles. 
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FIG. 4: Entanglement for distinguishable particles and iden- 
tical bosons versus to/t/o for = 4. Solid curve: E{p) for 
distinguishable particles. Dashed curve: Eb{p) for bosons. 

For large A^, it can be shown that n(x5) oc rik with X5 = 
hkt/m after a releasing time t. The correlation functions 
can then be derived analytically at to <C /7o and to ^ Uq 
using the wavefunctions ([2j) and (|3]). We find that 

C2{xb)\^gj, (X {ukTik') = h,Qh',QN{N - 1) (9) 

at to^Uo- The (5-function at /c = produces the major 
interference peaks corresponding to x = hGt/m with G = 
0, ±27r/ao in Fig.[2]and Fig.O At to <C Uo, for (fictitious) 
individual Mott insulator states, we have {ukTik') = {N — 
l)/N, which is a flat distribution over the /c-space; for 
the symmetric superposition of Mott insulator states, we 
have 

C2{xb)\i;sr {nknk^)^{N - 2 + 4,^/^'), (10) 

agreeing with that of (indistinguishable) bosons. Simi- 
lar results can be calculated for higher order correlation 
functions. For example, we have {nkUk'Tik") = 1/4 + 
'^^k,k'^k' ,k" for the symmetric state; while {nkUk'Tik") = 
3/8 for individual Mott insulator states. Therefore, the 
interference can be observed in the correlation functions 
at large N with reduced visibility. 

Finally, we study the entanglement in the ground 
states of both distinguishable particles and identical 
bosons. The statistics of the particles affects the en- 
tanglement in a very different way. For distinguishable 
particles, the entanglement can be characterized by cal- 
culating the entropy E{p) of a subsystem containing half 
the number of particles after tracing off the other half. 
The superfluid-like state at the strong tunneling limit is 
a separable state with E{p) =0. Starting from this state 
and gradually reducing the tunneling to zero, we obtain 
the symmetric superposition of the Mott insulator states 
which is a maximally entangled state with E{p) ^ 2.5, 
as is shown in Fig. [H Here, entanglement is generated 
because the interaction does not commute with the single 
particle Hamiltonian. 



For identical bosons, the study of entanglement is com- 
plicated by the fact that the particles are identical and 
the state space is much smaller than that for distinguish- 
able particles. We use an operational definition for en- 
tanglement presented in Ref. |^]. The entanglement can 
be described as the correlation between two subsystems 
that can be manipulated by local operations, with the 
definition 

N 

E,{p) = Y,PnE(Y[\4>g){^g\l[) (11) 

n=0 n n 

where Yin Pi"ojects the state to a subspace with n parti- 
cles in the left two sites. In Fig. [H we plot the entangle- 
ment for bosons. Instead of a monotonic increase with 
the interaction [/q as for the distinguishable particles, the 
entanglement for identical bosons £^5 reaches a peak at 
to/Uo ^ 0.16 and decreases to zero at to <C Uo where 
the projected state in the subspace is a pure state with 
no correlation with the state in the other subsystem and 
Ei, = E{c\cl\0){0\ciC2)=0. 

To conclude, we studied the ground state properties of 
the Hubbard model for distinguishable particles by exam- 
ining correlation functions and entanglement properties. 
In the strong interaction limit to <C /7o, the symmet- 
ric superposition of the Mott insulator states becomes 
the ground state and hence the quantum interference 
in the correlation functions and entanglement between 
the particles is found. Our results suggest that a bet- 
ter understanding of the role of quantum statistics in 
quantum phase transitions can be gained by studying 
distinguishable particles. Using a relatively small sys- 
tem of cold atoms [5, trapped ions[l2^, and exciton 
polariton[13, 14], theoretical predictions presented in this 
paper should be experimentally confirmed. 
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